Plot_Chat <- function (dat=msm1) {
  attach(dat)
  plot(1979:2010,obs_catch_1,ylim=c(0,1.5e6),ylab="Catch (t)",
       xlab="Year")
  lines(1979:2010,obs_catch_hat_1)
  points(1979:2010,obs_catch_2,ylim=c(0,1.5e6),pch=2)
  lines(1979:2010,obs_catch_hat_2)
  points(1979:2010,obs_catch_3,pch=3)
  lines(1979:2010,obs_catch_hat_3)
  detach(dat)
}
Plot_Chat(msm3)
  Plot_Bhat <- function (dat=msm1) {
    attach(dat)
    plot(yrs_srv_biom_1, srv_bio_hat_1,ylim=c(0,1e7),type="l",ylab="Pollock biomass",xlab="year")
    points(yrs_srv_biom_1, srv_bio_1,pch=19)
    plot(yrs_srv_biom_2, srv_bio_hat_2,ylim=c(0,1.5e6),type="l",ylab="Pacific cod biomass",xlab="year")
    points(yrs_srv_biom_2, srv_bio_2,pch=19)
    plot(yrs_srv_biom_3, srv_bio_hat_3,ylim=c(0,1e6),type="l",ylab="Arrowtooth biomass",xlab="year")
    points(yrs_srv_biom_3, srv_bio_3,pch=19)
    detach(dat)
  }
Plot_Bhat(msm3)